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We discuss nonspreading wave packets in one dimensional Schrodinger equation. We derive general 
rules for constructing nonspreading wave packets from a general potential V (#,£). The essential 
ingredients of a nonspreading wave packet, the shape function /(x), the motion d(t), the phase 
function (fi(x, t) are derived. Since the form of the shape of a nonspreading wave packet does not 
change in time, the shape equation should be time independent. We show that the shape function 
f(x) is the eigenfunction of the time independent Schrodinger equation with an effective potential 
V e fi and an energy E c s. We derive nonspreading wave packets found by Schrodinger, Senitzky, and 
Berry and Balazs as examples. We show that most stationary potentials can only support stationary 
nonspreading wave packets. We show how to construct moving nonspreading wave packets from time 
dependent potentials, which drive nonspreading wave packets into an arbitrary motion. 



PACS numbers: 03.65.-w, 03.65. Ge, 
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INTRODUCTION 



Nonspreading wave packets are those packets whose form of the probability density function ^(xji)! 2 under time 
evolution do not change. The shapes of energy eigenstates of a time independent potential are stationary and are 
therefore nonspreading in the trivial case. The first interesting nonspreading wave packets was found by Schrodinger 
[]]]. This packet is normalizabe, and is the shifted ground state of a harmonic oscillator. The motion of the packet 
behaves just like a classical particle shifted from its equilibrium position, executing the simple harmonic motion. After 
that, Senitzky extended the result showing that shifted higher energy eigenstates are also nonspreading Senitzky 
also showed that the expectation value of energy of a nonspreading wave packet is E n + E c i , where E n is the quantum 
mechanical energy of a wave packet when it is stationary, and E c i the classical energy of the particle. Thus, though 
the motion of a classical particle and a nonspreading wave packet is the same, their energy contents are different. We 
may view E n the structure energy of the stationary wave packet besides the classical energy. 

The extension of nonspreading packets from shifted ground state to shifted higher energy eigenstates were also 
found by Roy and Singh [3| , Yan 0] , and Mentrup et al. [B| , etc. 

Much later than Schrodinger, Berry and Balazs found another interesting nonspreading Airy packet Ai[x] in free 
space 0. Airy packets Ai[x] though bounded are not square integrable. An Airy packet does not move in uniform 
velocity in free space, instead, it accelerates. Thus, in contrast to Schrodinger packets, the motion of an Airy wave 
packet is different from that of a classical particle in free space. This is due to that Airy packet is not normalizable, 
and therefore does not really describe a particle. Nonspreading Airy packets even exist in time-dependent uniform 
force It had been shown that no other potentials except those described above will support nonspreading Airy 
packets Q. There are beginning a lot of investigation on the interesting optical Airy beams. The Airy optical beams 
was recently first observed by Christodoulides el. |8[ . Nonspreading wave packets in an imaginary potential was also 
observed by Stiizle el. ((J. These packets are called Michelangelo wave packets, which though reach a stationary 
width, the probability density however decays in time, and hence are not really form invariant. In second section, we 
show that it needs a real potential to have a form invariant packet. 

Nonspreading wave packets is an interesting and important subject. The essential ingredients of a nonspreading wave 
packet are its shape function f(x), its motion d(t), and its phase function (j)(x,t). From the previous works of many 
people, we can have a general discussion on this subject. We derive the general rule for constructing nonspreading 
wave packets in the second section. We derive Berry and Balazs's result in free space and a linear potential in the 
third section. We derive Schrodinger and Senitzky's result in the fourth section. We discuss results in other potentials 
in the fifth section showing that most stationary potentials can only support stationary wave packets. We discuss 
nonspreading wave packets in time dependent potentials in the final section. 



GENERAL DERIVATION 



We start with Schrodinger equation, 



dV(x,t) 
dt 



HV(x,t). 



(1) 



The Hamiltonian H is 



H = -^-d x d x + V{x,t). 



(2) 



The potential V(x,t) is allowed to be complex and is written as 



V(x,t) = V r {x,t) +i Vi(x,t), 



(3) 
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where V r (x, t) and Vi[x, t) are two real functions. We are looking for non-spreading wave packets in a general potential 
V(x, t). Let the initial wave packet be 

*(*,()) = /(as) e <9 <">, (4) 



where f(x) and 9(x) are real functions. We require the time evolved wave function ty(x,t) being of the form 

*(x,t) = f{q) e*<-**\ (5) 

where 

q = x-d(t). (6) 

The two functions d(t) and <p(x,t) are real, and satisfy the boundary conditions 
d(0) = 0. 

cf>(x,0)=6(x). (7) 



The function d(t) describes the motion of the packet, and d(t) is its group velocity. We have |^(x, 0)| 2 = f(x) 2 and 
|\&(a;,t)| 2 = /(g) 2 , so \^(x, t)\ 2 has the same form as \9{x, 0)| 2 ; therefore, ^(x,t) is a nonspreading wave packet. 

We need to determine the functional form of f(q) and the related phase function <f>{x, t). We rewrite the Schrodinger 
equation as 

d^lx t) h 2 

Substituting |5} into {SJ and separating the equation into real and imaginary parts, we have then two equations. The 
real part gives the equation: 

|-/"(<?) = [ -V r (x,t) + W t ct>{x,t) + h2 ^ x ^ 2 ]f{q) . (9) 
Im Ira 



The imaginary part gives the equation: 

+ —d x <f>(x,t)] /'(g) = [ Vi(x,t) - ^-d x d x cj>(x,t) } /(g). (10) 

We have used the notation: /'(g) represents the derivative of /(g) with respect to g, and d(t) the derivative of d(t) 
with respect to t. The real part equation relates /"(g) and /(g). The imaginary part equation relates /'(g) and /(g). 
Eq. (|10p . when multiplied by 2m/ (g), can be rearranged into a form as a flux equation 



d x j(x,t) = 2mV l (x,t)f(q) 2 , 



(11) 



4 



where 

j(x, t) = [-mh d(t) + h 2 d x <f>(x, t)} f{q) 2 . (12) 
Taking integration of x of both sides of (jllp from — oo to oo, we have 

POO 

j(oo,t) - i(-oo,*) =/ 2mV l (x,t)f(q) 2 dx. (13) 



We consider bounded f(x), such that f(x) — > for x — > oo. Hence the left hand side is zero. The right hand 
side, if Vi is not a zero function, in general is not zero. This means a complex potential in general does not support 
nonspreading packets. It is then to require that 

Vi(x,t)=0. (14) 



From (fTTjl . we have 

8 x j(x,t) = 0. (15) 



(fT5|l shows that j(x, t) can only be a function of time. Then 

j(x,t) - [-mhd{t) + h 2 d x <t>{x,t)] f 2 (q) = b{t). (16) 

Using again the boundary condition of f(x) at infinity, we have 

b(t) = 0. (17) 

From (|T6|) and (fl7j) . we then have the important result 

-mhd(t) + h 2 d x (j)(x, t) = 0. (18) 

From (fl~8|) . the form of <j>(x,t) is determined to be 

4>(x,t)=<l> 1 (t)x + <t>o(t). (19) 
Mt) = (20) 
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Thus the phase function, <f>(x,t), is restricted to be a function linear in x. From (114[) . the potential is real, we then 
simply write the potential as V(x,t). 

We next need to determine the functional form of f(x). Substituting (flT))) . (f2T))) into ©, we have the equation 
determining the shape of the function f(q). 

ti 2 md(t) 2 

— f"(q) = [ V(x, t) + — i^- + md(t) x + h<p (t) ]/(<?)• (21) 
Zm Z 



The left-hand side of ([2"T]) is a function of q, so should be the right-hand side. We replace the variable x by q + d(t), 
and then expand the terms in the bracket of the right hand side in terms of the powers of q. We have 

fe2 00 

f'\ q ) = C£C n q n )f{q). (22) 
Zm e - — ' 

n=0 



Intuitively, each of the coefficient C n is time dependent constant, and should be written as C n (t). However, in order 
to be consistent with the left-hand side being a function of q only, C n can only be time independent constant. We call 
the requirement that C n should be time independent, the consistency condition. The constant Co contains terms 
md(t) 2 / 2 + md(t) d(t)+h<po(t), and also the constant term from V(x, t), which is, V(d(t),t). We denote the coefficient 
Co by — -E c ff and denote all the rest terms in the summation by V e g. Then 

-E eS = V(d(t),t) + + md(t)d(t) + HMt)- (23) 



and 

V cS {q) = V(q + d(t),t) - V(d(t), t) + md\t) q. (24) 

Eq. (|2ip can then be rewritten as a form similar to that of the time independent Schrodinger equation. That is 

-^f"(Q) + V eS (q)f(q) =E eS f(q). (25) 



We call (|25|) the shape equation. We see that the function f(q) is the eigenfunction of a Hamiltonian with a potential 
V c ff(q). Since the functional form of a nonspreading wave packet does not change in time, the shape equation being 
related to the time independent Schrodinger equation is quite interesting. Finally, by setting E e g to one of the energy 
eigenvalues, we can solve 4>o(t) from (|23|) . 

We have obtained the general rule for constructing nonspreading wave packets from a general potential V(x, t). We 
list the following formulas for a conclusion: 
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9(x,0) =f(x) e l6(x \ (26) 

9(x,t) = /(g) e^l (27) 

q = x-d(t). (28) 

cf>(x,t) =</>i(t) x + fo(t)- (29) 

= (30) 

" 2^ r(<z) + Kff(9) = EcS f{q) - (31) 

Kfr(g) = V(g + d(t), *) - V{d(t), t) + md{t) g. (32) 

E cS = -[ V(d(t), t) + + d W + h Mt) ], (33) 



together with the consistency condition, and the boundary condition: /(a;) — > for |x| — > oc. 
In what follows, we give applications to these formulas. 

BERRY AND BALAZS'S RESULTS 

We start from a real potential V(x, t) — 0. To determine the shapes of the nonspreading wave packets in this free 
space, we first calculate the effective potential V c g(q). We have from 



V cS (q) = md{t) q. (34) 

Since V e s can only be a function of q, it needs that d(t) be a constant. We set 

d\t) = B 3 /{2m 2 ). (35) 

Then V c g — B 3 q/(2m), which is a linear potential. Substituting this into pip , we obtain the shape equation 

/"(g) - ^ g /(g) = =?^m (36) 



The solution of (|36[) is the Airy function. Thus the nonspreading wave packet in free space is determined to be Airy 
functions. We consider here the simplest solution by taking E e s = 0. The solution of Eq. (f3"6]) is then Ai[(B/h 2 ^ 3 ) q], 
where we neglect the divergent solution Bi[(B /h 2 / 3 ) q]. 

The motion of the packet, d(t), can be solved from ([33)) . Choosing the boundary conditions as d(0) — 0, and 
d(0) = 0, we have the simplest solution: 



(37) 
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We obtain ^(t) from ([30]). Then 

Mt) = |£. (38) 

We obtain 4>o{t) from (|3"3"|) with .E eff =0, and the boundary condition 4>o(0) = 0, we have 



We compare these results with Berry's result: 



We see they are the same. From (|36|) . we also see that the only nonspreading wave packets in free space is the Airy 
function [l(|. The trivial plane wave solution corresponds to taking B = in (|36|) . 
We next consider the real potential 

V(x,t) = -f(t)a;. (41) 

From (f3"2")h we have 

y off ((?) = [-F(t) + md(t)] 9- (42) 

For consistency, it needs that \—F(f) + md(t)) be a constant. We set 

-F(t) + md\t) = B 3 /(2m). (43) 



Then V c s = B 3 /(2m)q. Thus the nonspreading wave packets in the potential V(x, t) — —F(t)x are the Airy functions 
again, as found by Berry and Balazs. Again, we calculate 
The rest of the complete results are easily shown to be 

BH 2 1 '* r T 
lo Jo 



d(t) = — j + — / / F(x)dxdT. (44) 



m = ^{ + \ f F{r)dr. (45) 
2mn h J Q 



ft i<t rt i-t ?x 



2m 2 h 



pt pT pt pT pX 

[ I t I F(x)dxdr+ / / / F(y)dydxdr}. 
Jo Jo Jo Jo Jo 



(46) 
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Using the identity 

[ {t- T)F{r)dT = [ [ F{x)dxdT. (47) 

Jo JO JQ 

we see these results are the same as those of Berry and Balazs's @. 

SCHRODINGER AND SENITZKY'S RESULTS 

We consider the potential 

V(x) = ^muj 2 x 2 . (48) 

From (|32|) . we have 

V eS (q) = ]^nu 2 q 2 + m[ d\t) + uj 2 d(t) } q. (49) 

It needs that d + cu 2 d(t) be a constant. For a simplest solution, we set the constant to be zero. Hence 

d+ uj 2 d(t) = 0. (50) 



Then d(t) executes simple harmonic motion. This is the same as the classical solution of x(t). Taking (f50|) into (|49|) . 
we have 

Vctf(q) = X -mu?q 2 . (51) 



The shape equation then is of the form 

~f"(q) + \™^ 2 q 2 f(q) = E cS f(q) (52) 
2m 2 



We see f(q) = t() n (q), the energy eigenfunctions of the simple harmonic oscillator. Each eigenfunction ip n (q) then 
offers a nonspreading wave packet. These packets though with the shape of stationary energy eigenfunction, however, 
they are moving and are in fact executing the simple harmonic motion. We have then reproduced Schrodinger and 
Senitzky 's result. 

Starting from (fB"0"|) , we have the following results 



d(t) — xq sin(uit). 



(53) 
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h x {t) = cos(ut). (54) 



6 (t) = ^ sin(2wt) - -j-. (55) 



We also easily see that if the frequency is time dependent, that is u = ui(t), then the effective potential from (|19")) is 
Vcs(q) = ^mu){tfq 2 + m{ d(t) + uj(t) 2 d(t) } q. (56) 

From the consistency requirement, we need the coefficients of q and q 2 be constant, hence ui should be a constant. 
Thus a harmonic oscillator with a time dependent frequency can not support a nonspreading wave packet, as described 
in Yan's paper 

POTENTIAL V(x) = \x n ,n > 3 

For potentials of the form 

V(x) = Xx n ,n > 3, (57) 

it is trivial to see that from V c s and the consistency equation, the motion d{t) must be a constant, function 
From (|32|) . we have 

V cS (q) = \[q + d(t)] n - Xd(t) n + md q. (58) 

Expanding X[q + d(t)] n in terms of powers of q, we have the first few terms 

Xq n + Xud^q 11 - 1 + ■■■ . (59) 



The second term of (|59p shows that d(t) needs to be a constant. Thus the packet does not move. This implies that 
the nonspreading packets can only be stationary packets. We set this constant as d, then we have 

d(t) = d. (60) 



And then we have 



V cS {q) = X[q + d)] n -Xd n . 



(61) 
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Substituting this to (|3"Tj) , we have the shape equation of the form 

~f"(q) + % + d)] n f(q) = Ef(q), (62) 
Zm 

where E = E e ft + Xd n . As q + d = x, formula (|6"2")l shows that the function f{q) = ip n (x), the energy eigenfunction 
of the Hamiltonian with potential V(x) = Ax™. Since the packet is the energy eigenfunction, it is a stationary wave; 
and therefore it does not move. We conclude that there is no moving nonspreading packet in a potential of the form 
V(x,t) — Xx n ,n > 3. The same argument applies also to showing that a potential, such as V(x) = X/x, does not 
support moving nonspreading wave packets. We see most stationary potentials can only support stationary wave 
packets. 



TIME DEPENDENT POTENTIALS WITH MOVING NONSPREADING WAVE PACKETS 

From above results, hence, in order to have moving nonspreading wave packets, we need to consider time dependent 
potentials. The simplest way is to move the potential. That is to consider V(x — d(t)). In this way, there are in fact 
many nonspreading wave packets with arbitrary motion that can be constructed. For instance, we start from a given 
motion, such as d(t) = < 2 , etc., and then consider the time dependent potential 

V(x,t) = —mui 2 (x — d(t)) 2 — mxd(t) — — muj 2 x 2 sin(ojt) 2 . (63) 

The effective potential calculated from (f6"3"| and (|52"]) is 

V cS (q) = \mu?q 2 . (64) 



Thus, the potential in Eq. (|63|) offers nonspreading wave packets, with the shapes of the eigenfunctions of SHO, 
moving with the arbitrary function d(t). If we set d(t) = xq sin{u)t), and substituting this to (|63p . we find the needed 
potential for this motion is V(x,t) = ^muj 2 x 2 , as expected. 
Another example is considering the time dependent potential 

V(x, t) = X(x - d(t)) 4 - mxd\t). (65) 

The effective potential can be calculated to be 

V cS (q) = X q\ (66) 



The potential in Eq. (|6"5|) then offers nonspreading wave packets with the shapes of the eigenfunctions of the potential 
V(x) — Ax 4 . These nonspreading wave packets move with the arbitrary function d(t) as the moving potential does. 
We see that in order to move a nonspreading wave packet, we need a moving potential, and also a linear potential 
—mxd(t), which offers a force md(t); this is just the force needed for a particle with a motion d(t). In conclusion, a 
moving potential A(x — c?(t)) 4 , together with a time dependent linear potential —mxd(t), will drive stationary-shaped 
packets going around with a journey d(t). It is like to have water moving in space, we put it in a bowl, and carry the 
bowl going around. 
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